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Abstract 

We use our Nother Symmetry Approach to study the Einstein equations mini- 
mally coupled with a scalar field, in the case of Bianchi universes of class A and B. 
Possible cases, when such symmetries exist, are found and two examples of exact 
integration of the equations of motion are given in the cases of Bianchi AI and BV. 
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1 Introduction 



Astronomical observations imply that, on a very large scale, our universe can be simply 
considered as homogeneous and isotropic for cosmological purposes, being thus well de- 
scribed by a Friedman-Robertson- Walker (FRW) metric IJ. It was suggested @ that 
such a scenario could be nothing else but the result of a different previous history and, 
therefore, be also consistent with a large amount of primordial anisotropy. 

If, moreover, one adopts an inflationary scenario to describe the very early history of 
the universe, it can be expected that what is actually observed today does not depend 
on specific initial conditions, including possible initial anisotropics 0. On the contrary, 
pre-inflationary anisotropies can play a role if primordial friction coefficients, responsible 
for the coupling between the gravitational field and the inflaton field, are supposed. As 
a matter of fact, such coefficients are proportional to the rate of expansion of volume, 
and their increasing makes inflation work better (see for example [|J). 

Furthermore, since we generally expect that all homogeneous anisotropies do not in- 
crease, it is not strange that we find very small anisotropies in the microwave background 
(10~' -T- 10~ 4 ) on large and small angular scales [0, 0, ||. Anyway, it is possible to con- 
jecture that anisotropic cosmological models isotropize as time approaches our epoch ||, 
0, which makes it meaningful to study situations in which isotropy is not necessary at 
early times. 

On the other side, as we just mentioned before, introduction of a scalar field in 
cosmology is widely used, to obtain inflationary eras at the very beginning and to treat 
primordial matter fields in the schemes of Grand Unified Theories. 

For all these reasons, it is worthwhile to understand dynamics of anisotropy and to 
search for anisotropic cosmological models which are exactly solvable. 

In this paper, we investigate homogeneous theories of gravity, in which a scalar field 
is minimally coupled to gravity, searching for point symmetries in the cosmological La- 
grangian density which allow to solve exactly the dynamical problem (we begun such a 
program in M). 

We start from the action 

A = J d'x^-g (r + \g^p-^ - %)) , (1) 

where R is the Ricci scalar curvature for a spatially homogeneous spacetime, described 
by the metric 

ds 2 = -dt 2 + g ab (t)uj a LU b , (2) 

u a being 1-forms, not necessarily closed (a,b = 1, 2, 3). By (/i = 0, 1, 2, 3) we denote 
the covariant derivative of the scalar field ip; g is the determinant of the metric, and V(tp) 
is the potential for <p. 

Brans-Dicke gravitational theories have been recently M examined in a wider context. 
Anyway, a search for point symmetries in vacuum situations, in which <p is nonminimally 
coupled to gravity as well, has been recently performed in [|Kj in a form, however, which 
is not immediately reducible to the minimally coupled ones. 
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Here we search for Nother symmetries of Lagrangian ([I]), in absence of ordinary 
matter, in order to select all Bianchi models which are exactly solvable by the so called 



Nother Symmetry Approach |8j, [I I 



The paper is organized as follows. The Bianchi classification and the Lagrangian 
description are sketched in Sect. 2. In Sect. 3 we apply our procedure to determine the 
cases where symmetries are present. Sect. 4 is devoted to the analysis of the results, and in 
Sect. 5 we give two examples of exact integration of the equations of motion. Conclusions 
are drawn in Sect. 6. 

2 Bianchi classification and Lagrangian formulation 

Bianchi spacetimes, in modern geometric language, are G 3 — principal fiber bundles over 
R. The bundle being trivializable is usually written as M 4 = R x G 3 . G 3 stays for 
the three-dimensional group of isometries of the metric on M 4 . If one is not paying 
attention to the topology of G3 (i.e. not distinguishing between the simply connected 
group with respect to the quotient groups obtained by quotienting it with discrete normal 
subgroups (not distinguishing among Lie groups with the same Lie algebra) we can use a 
parallelization of M 4 with the vector fields X , X%, X 2 , X 3 with X the "time" -generator 
and X\, X 2 , X 3 the Killing vector fields of the metric closing on the Lie algebra of G3. Of 
course, when X l5 X 2 ,X 3 are known (and complete) we can integrate them to the action 
of G 3 on M 4 . 

In the framework of this procedure, in 1897 Bianchi worked out the first classifica- 
tion of three-dimensional Lie algebras |I3|. Of course, metrics for which corresponding 



groups are isometries describe cosmological models with specified symmetry properties. 
Bianchi started with the adjoint representation and his approach has been followed since 
then. Here, we would like to comment briefly on the classification by using the coadjoint 
representation, as recently proposed in |14| and |15 . 



As it is well known a Lie group is locally characterized by its structure constants c bc . 
If (X a ) is a basis for the vector space underlying the algebra Q we have [X a , X b ] = c c ab X c 
and the Jacobi identity [X a , [X b , X c }} + [X b , [X c , X a }} + [X c , [X a , X b }} = gives a quadratic 
constraint on the structure constants. 

By identifying the vector space Q with the linear functions on Q*, Q = Lin((?*,R), 
we can define a Poisson bracket on Q* by setting {x a , x b } = c c ab x c . The classification 
problem is reduced now to the problem of classifying linear Poisson brackets on a three- 
dimensional vector space. By using the transformation properties of a Poisson bracket 
under change of coordinates one shows that with any bracket is associated a bivector 
field A. In our case A = c c ab x c -^- A 

In three dimensions, any bivector field is associated with a one-form a by means of 
the contraction with a volume form Q, say Q = dx\/\dx2/\dx 3 . Thus, we have to consider 
and classify quadratic one-forms a = c c ab x c e abr dx r = H cr x c dx r . The Jacobi identity is 



equivalent to a A da = [14 



The one-form a can be written by decomposing H cr into symmetric and skew-symmetric 
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part 

a = -H^ cr \x r dx c + x c dx r ) + -H^ cr \x r dx c — x c dx r ) . (3) 

By using linear transformations on the three dimensional vector space, we can bring a 
to the normal form 

a = -d(Nix\ + N 2 x\ + A^Xg) + ^(x 2 dx 3 - x 3 dx 2 ) ■ (4) 

Imposing the Jacobi identity, we find aNi = and, by means of a suitable change of 
coordinate basis, we may set Ni = 0, —1, +1. Thus all quadratic one-forms giving rise to 
three dimensional Lie algebras are classified in two classes 

1. Bianchi A , a exact, i.e. a = 

2. Bianchi B, da ^ 0, i.e. a / 0. 

In this approach, class A is associated with a foliation of R 3 given by the Casimir 
function N\x\ + N 2 x 2 + N 3 x\. The level sets of this function are the coadjoint orbits 
of the Lie group G 3 associated with the Lie algebra Q . They give information on the 
topology of the group. Clearly, if N\N 2 N 3 ^ 0, we have either 0(3) or 0(2, 1), according 
to the signs of NiNj (the second, if one of them is negative). If at least one of the A^ is 
zero, the group is noncompact. 

The parameters a, A, are sufficient to characterize Bianchi classification, according to 
the following table 

Table I-Classification of Bianchi types 



Type 


a 


Nt 


N 2 


N 3 


I 














II 





1 








III 


1 





1 


-1 


IV 


1 








1 


V 


1 











VI 





1 


-1 





VI,, 


a 





1 


-1 


VII 





1 


1 





Vllfc 


a 





1 


1 


VIII 





1 


1 


-1 


IX 





1 


1 


1 



A parameter h = a 2 /(N 2 N 3 ), is also introduced for sub classifying tipes VI and VII. 

Because any Lie group is parallelizable, we can use left-invariant (or right invariant) 
one-forms in the metric (||D, so that it will be invariant under left (or right) action of G3 
on M 4 . By using a diagonal form for g a b and selecting (for obvious reasons) a positive 
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definite one, we can parametrize g a b by factoring det ||<7 &||, namely, following Misner, we 
obtain [0, 0, 

g ab (t)=e 2 ^e 2 ^ , (5) 

where the matrix (5 is diagonal and traceless, thus depending on two variables only. A 
widely used choice is |L2[ 



diag f-ift + ^(3 2 , -\fS x - ^(3 2 , fa 



(6) 



Being y/—g = e , the expansion of volume is wholly determined by A. As to the 
shear, it is determined only by the This is typical only for the 'orthogonal' case, 
when there are no effects of rotation and tilt. Spatially homogeneous sections are then 
orthogonal to the fluid flow vector in the universe [|12j ||16|| . 



As shown in |12j , the second-order Einstein equations of all Bianchi types, except IV 
and VII/^o, can be derived by a point Lagrangian. In Sect. 5, we shall consider also the 
first order Einstein equation, which is always a first integral of the second order ones. 

Then the general expression for the Lagrangian is the following 



C 



.,3A 



R * + 6 >?--(fi + pl)-^ + 2V(tp) 



(7) 



with the scalar curvature of the spacelike hypersurfaces given by 



R* 



-2A 



for the types belonging to class A, and by 

R* = 2aV 2A f °> 

with 



+ -N^N^l + NiN 2 N 3 ) 



N * N * \ -0 



(Ar 2 AT 3 /5i + J-3a?N 2 N 3 p 2 



(9) 



(10) 



3a 2 - N 2 N 3 

for types in class B, with n a a = 0. 

As said above, for types IV and VII^o> the Lagrangian is not known, and therefore 
we shall not consider them here. 



Nother symmetries 



According to a well established procedure |]q|, |[Llj| , we look for Nother symmetries of the 



Lagrangian (|7|). In other words we look for a vector field X on the configuration space, 
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such that the Lie derivative with respect to the tangent lift X T of X on the tangent space 
vanishes 

L x tC = . (11) 

The conditions imposed by this equation generally puts restrictions on the class of allowed 
potentials for ip. 

But let us first observe that there are a number of trivial cases. 



• If V = const., C does not depend on tp, which is therefore cyclic. Thus, for all 
classes and types, we always have the symmetry 

x - £ ■ < 12 > 

• In Bianchi AI and BV, we have that R* does not depend on /5 2 , which are 
therefore cyclic. We have thus the symmetries 

• x ^k- (13) 

Moreover, since the term f}\ + /3| is rotationally invariant, we have the symmetry 

*= A ilr A 4 +A 4- A £- (14) 

• In Bianchi All, R* depends on (3\ only, so that we have the symmetry 

x= m- (15) 

In the following, using (|Tl|), we investigate the possibility of other symmetries, not so 
immediately evident. 

First of all, we have to say that the configuration space is made of four variables, 
namely A, Pi, /3 2 and <p. Let us set 

d d d d 

x - L lTx +B ^ + B ^ + F T^ < 16 » 

where L, B\, B2 and F are unknown functions of A, (3\, (32, <p. The tangent lift will be 
_ T v dL d dB 1 d dB 2 d dF d 

j^i _ j£ j 1 i 1 f. 1 (17) 

dt d\ dt d$i dt d$ 2 dt dp 1 
where d/dt means Lie derivative along the dynamical vector field; for example, it is 

dL dL • dL ■ dL ■ dL 

Ht = a\ + W/ 1 + + ^ ' (18) 



The expression L x tC is a homogeneous quadratic polynomial in A, Pi, p 2 , (p, plus a 
term of zeroth degree. Therefore, it vanishes if and only if all the coefficients are zero 
independently. This leads to a system of eleven partial differential equations 



dL 

OA 


= ; 


dBi 

3L + 2— - = 

dpi 


(19) 


dB 2 
dp 2 


= 


dF 

dip 


(20) 


dL dB l 
dPi dX 


= 


dL dB 2 _ 
dp 2 dX 


(21) 


dip dX 


= ; 


dB 1 dB 2 
dp 2 dp x 


(22) 


dB 1 dF 

6 dip + dp, 


= 


dB 2 dF 
3— ^ + 2— = 

dtp dp 2 


(23) 



L (R* + 6V(<p) + NiN 2 N 3 (l + N^Ns)) + B 1 — + B 2 —- + 2FV' = , (24) 

op\ op 2 

which are valid for both classes A and B, with relatively different expressions for R* . 

We first obtain a general solution for equations ([19| - p3D , and then use (0) as a con- 
straint to limit the class of solutions and find a condition on the potential. 

From Eq. (|19]a) , we immediately get 



so that from (pi 



and from (|22|a) 



B±o = — 



L = e-l x L((3 u (3 2 ,ip) , 
8 dL 



--—e-^ + B^^fo^) 

6 Opi o 



dL 

F = -4^e--i x + F(p 1 ,(3 2 , V ) 



Using fl2"2]b), we get then 

_ 8 d 2 L 



dip 



3 A dBi 8 d 2 L 



dp 2 3d(hdp 2 e 2 
Since Bi t2 do not depend on A, it must be 



dB 2 
dPx 



(25) 
(26) 
(27) 

(28) 



so that 

L = L 1 {(3 l ^) + L 2 {(3 2lV ) ■ (30) 

moreover fl22| b) gives 

dB 1 dB 2 



d(3 2 ' df3 1 

Substituting into (|T9|b), we get 



,«,/= - 16d 2 LA dB 1 

3e - 3/x [ L i +L >-^W + wr°- 



(31) 



(32) 



Now, since only L 2 and B 1 depend on /3 2 , and B 1 does not depend on A, it must be 

OB 

1^0 ; ^ = 0. (33) 
Substituting into (|20|a) we get analogously 

L, = ; ^ = 0. (34) 



Thus, we find that L = and i?i = Bi((3 2 ,tp), B 2 = B 2 (f3i,(p) . From (|22"|b) we obtain 



now 



B x = f{<p)fo + 9l (<p) ; B 2 = -f(<p)pi + g 2 (<fi) , (35) 
while from (pOfc) and (7p2^a) we get F = F(j3i,j3 2 ), and from (f23|a) 

3/'/3 2 + 3g[ = -2— . (36) 

The left hand side does not depend on fli, thus 

F = -ltf 1 -f%)Pi + F(p 2 ) , (37) 



and, analogously, from 



F = -l(g' 2 -f(3 1 )(3 2 + F((3 1 ) . (38) 



By comparison we get 

f = ; g" = g 2 =0, (39) 

- ~~ F(3 X = -F ; ^ 2 + F/5 2 = F , (40) 

where F is an arbitrary constant. The general solution is thus 

L = ; B x = c/3 2 + c x ip + c (41) 

5 2 = -cPi + c 2 y + c' ; F — F - ~(ciPi + c 2 /3 2 ) , (42) 



with c, ci, c 2 , Co, c' , F arbitrary constants. 

We have now to plug solutions (0) into ([P]) and check for compatibility. Being L = 

we get 

B € + *if = -^- <«» 

It is easy to see that the two sides must be zero independently, so that we have two 
subcases: V = const, and F — 0. In both of them (pE3|) is satisfied for types I and V. 
Examining the other types, we find that for Class A Eq . fl43|) is verified only for type II 
and B\ = 0, while for Class B Eq. ([43|) is verified for type Vlh^o iff 

1.1 

c = ; c 2 = /=ci ; c = y=c . (44) 

ay 6 ay 6 

The other types are excluded from the beginning, as said above. 

4 Results 

Case 1 — V = costant 

The situation for V = const, is summarized as follows 
la) Bianchi AI 

The most general symmetry is written (on Q) as 

d d 3 d 

X = (cfo + Cl tp + c )— + (-c/3i + c 2 y + c' )— + (F - -(ciA + c 2 (3 2 )) — , (45) 

with obvious lift to TQ. A basis of symmetries on TQ is given by 

x i = ;Jr ; x * = t!r ; x ^ = ^ (46) 

5 3 ^ 5 • 5 3- <9 
X 4 = y?— - -A— + ^-^ - -A— 47 

53^9 <9 3 A <9 , . 

d(3 2 2 dy? df3 2 2 <9y? 

<9/3 2 96 df3 2 d/3i 

It turns out that only five of them are independent. For instance, we can take the first 
five of them. With these symmetries, we associate the following constants of the motion 

iri = e 3A A ; K 2 = e 3X fo ; K 3 = e 3 ^ (50) 

K, = e 3A (ip$i - ; K 5 = e 3A - ~(3 2 <p) ; K 6 = c 3 \[lj 2 - fafo) . 

(51) 



These symmetries close on the following Lie algebra 

[X i ,X j ]=0 , i,j = 1,2,3 ; [X 1 ,X 4 ] = ~X 3 ; [X 2 ,X 4 ]=0 (52) 

[X u X s ] = ; [X 2 ,X 5 ] = ~X 3 ; [X 3 ,X i ]=X 1 ; [X 3 ,X 5 ]=X 2 (53) 

[X 4 ,X 5 ] = ^X 6 ; [X lt X 6 ] = X 2 ; [X 2 ,X 6 ] = -X 1 (54) 

[X 3 ,X 6 }=0 ; [X 4 ,X 6 ]=X 5 ; [X 5 ,X 6 ]=X 4 , (55) 

and the close on the same algebra in terms of Poisson brackets. 

This situation was examined in || (with different variables). In that paper, it was 
shown how to use the symmetries in order to obtain exact integration of the dynamics. 
We give an example below, 
lb) Bianchi All 

The general form of X (on Q) is 

*=(<W + *)^ + (*-f<>ft)£. «*> 
from which we get the independent fields (on TQ) 

9 v d v d 3 d . d 3- d 

X x = — ; X 2 = — ; X 3 = p— - -fa— + tp—r- - -p 2 — , 57 

df3 2 dp d[3 2 2 dip df3 2 2 dp 

with Lie algebra 

[X 1; X 2 ] = ; [X U X 3 ] = -^X 2 ; [X 2 , X 3 ] = X, . (58) 

The associated constants of the motion are a subset of the above ones. 

In spite of the existence of three symmetries, the exact integration is not straitforward 
in this case, so that we shall not study the situation here, 
lc) Bianchi BV 

The situation is quite similar to Bianchi AI , with the same symmetries and constants 
of the motion. The equations of course are different. An example of exact integration is 
given below. 
Id) Bianchi BVI^ 



We have 



- V| — : AV -.- • -Wttt- ( 59 ) 

aV3 o(3 2 

3 / , 3 , \ d 

X 3 = p— tev^t - o A-nfthr + 



d 


■ x - 9 






dp 




d 


1 d 




aV3^df3 2 


. d 


1 . d 


p> — — 


ay/3 d$ 2 



dp 
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with algebra 



[X 1: X 2 ] = ; [Xt, X 3 ) = X 2 ■ [X 2 , X 3 ] = ( -L _ 1) X 2 (61) 



3a 2 



and constants of the motion 



K 1 = e^<p ; K 2 = e 3A (ft - ^/? 2 ) (62) 

#3 = e 3A ^-3^i + ^/3 2 + 3/3^ - . (63) 

Case 2 — V 7^ const. 
2a) Bianchi AI 

We obtain the following independent symmetries 

d d d d ■ d ■ d 

X 1 = — ; X 2 = — ■ X 3 = p 1 —-p 2 —+p 1 - F P 2 - F (64) 
<9A 9#l df3 2 dpi 

with algebra 

[X 1? X 2 ]=0 ; [X 1 ,X 3 ] = -X 2 ; [X 2> X 3 ] = ^ (65) 
and constants of the motion 

K 1 = e 3X f3 1 ; K 2 = e 3X $ 2 ; K 3 = e 3A (/? 2( #i - /?i/3 2 ) . (66) 

This case was examined already in || and, although it was not possible to achieve 
exact integration, qualitative consideration on the behaviour of the solutions were made 
possible. 

2b) Bianchi All 

Only one simmetry is left, namely 

x =h- (67) 

2c) Bianchi BV 

Again the situation is essentially the same as for Bianchi AI. 
2d) Bianchi BVII 

We have only one symmetry 

x = m ~ m (68) 



with the corresponding constant of the motion 

"i-n/3 



K = e 3A (fa - -^=P 2 ) . (69) 
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5 Examples of exact integration 



In this Section, as an example, we obtain exact integration in the Bianchi AI and BV 
cases, when the potential V(cp) is zero. 

As said earlier, the first order Einstein equation is a first integral of the second order 
equations derived from the Lagrangian (0). It turns out that it is always equivalent to 

E c d ^^-q*-£ = 0. (70) 

It is in fact this equation that is left after reduction. Substituting the constants of the 
motion found above we get 

6A 2 - K 2 e~ 6X - a 2 e~ 2X = , (71) 

where K 2 = 3(K 2 + K 2 )/2 + K 2 and a 2 = Qa 2 e^. 

In Bianchi AI, i.e. when a = 0, we obtain immediate integration 



A = - lot 

3 f 



\K(t - to) 



(72) 



which coincides with the one found in || and therefore, we avoid further comments on 
it. In Bianchi BV, setting x = e 2X , we obtain the equation 

3 K 2 

±& _ £L _ o? = o , (73) 
2 x 

which integrates to 

[2 



g t = ^(1 +x) -log(Vz + VTTx) , (74) 

where we have set K — 1 (that is x(0) = 0), and a = 1 (that is, a rescaling of x and t). 

This relation cannot be inverted but can be used for a qualitative analysis. Indeed, 
it is easy to see that, for small t, and hence small x, we have that the mean scale factor 
e A is proportional to t 1 ^ 2 , while for large t it is proportional to t 1 / 3 . It is also possible to 
see that power law-inflation never occours. 

As for the behaviour of (p, it is immediately derivable from the expression of K3. In 
the first case we get 

<P = - to) > (75) 

while in the second we have ip oc t~ x l 2 for small t and <p oc logt for large t. 
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6 Conclusions 



We have found all Nother symmetries for the Bianchi universes for which a Lagrangian 
function, of the type given in is known. This allowed immediate exact integration 



in simple cases. In the other cases we have seen that the number of symmetries is often 
sufficiently high to permit a good reduction of the configuration space. These situation 
are of great physical interest and a complete analysis shall be done elsewhere. 

It is important to observe that the Lagrangian (0), in the case when we want to add 
ordinary matter (in form of dust decoupled from the scalar field), changes only by an 
additive constant. This is not irrelevant, because this constant changes the equation 
E c = 0, but of course the second order equations and the structure of the symmetries 
are exactly the same. This means that all the discussion above is valid for this case, 
except for the integration, which is now more complicate. 

Other types of ordinary matter only give a contribution in the coefficient of L in Eq. 
([24]). Since L = 0, we see that again the structure of symmetries is the same and the 
problems for integration may still arise from the energy condition. 
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